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On the inversion of y^eV in terms of associated Stirling numbers 

by 

D. J. Jeffrey\ R. M. Corless\ D. E. G. Hare^, and D. E. Knuth^ 

The function y = $a(a^), the solution of y°'e^ = x for x and y large enough, 
has a series expansion in terms of Inx and Inlnx, with coefficients given in terms 
of Stirling cycle numbers. It is shown that this expansion converges for x > (ae)" 
for a > 1. It is also shown that new expansions can be obtained for $q. in terms 
of associated Stirling numbers. The new expansions converge more rapidly and 
on a larger domain. 

1. Stirling numbers — Stirling cycle numbers [J^] are defined [1] by 



[la) + z) =m\J2i 
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The numbers (— l)"^"*"'" are also called Stirling numbers of the first kind [8]. Stirling 
subset numbers {^}, also called Stirling numbers of the second kind, are defined by 
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(16) (e^_l)- = ^!^jMi^ 



and 2-associated Stirling subset numbers {J^}>2 defined by [2, exercise 5.7; 7, p. 296; 
9, §4.5] 



;ic) (e^-l-^)^ = m!^ 
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2. Solution by Comtet of y^e^ = x. — The fixed real a, we let ^a{x) be the value 
of y that is the unique positive solution of the equation y"e^ = x. If a is negative, then 
y > —a and x > e~"(— a)". An asymptotic expansion for $q,(x), in terms of Stirling cycle 
numbers and the quantities Li = Inx and L2 — Inlnx, is given in the following theorem 
[3,5]. 

Theorem 1. — With the preceding notation, the function $q,(x) has the following 
series development, convergent if x is large enough. 



(2a) $«(x) = Li - aL2 + a V ^ V (-1^-+- 
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Proof. - We recall some details of the proof given in [5] for use below. We introduce 
a function w{x) defined by 

{2b) y = $a(a^) = Li — aL2 + aw , 
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which satisfies 

(2c) l-e-"" + aw-T = 0, (t=-^, t = = aln (-) . 

Li Li \(7/ 

By the Lagrange Inversion Theorem [4], w has the expansion 
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One converts from a and r back to Li and L2 to complete the theorem. 

Since the domain of convergence of (2a) is described only as large enough' by 
de Bruijn and Comtet, we give a stronger statement in the next theorem. 

Theorem 2. — For a > 1, the series (2a) is convergent for x > (ae)'^, while for a < 1 
it is convergent for x > e. 

Abbreviated Proof. - We let f{w) = aw — r and g{w) = 1 — e~'^ . For a > 1 and 
X > (ae)", we define 5>0by5 = l — ln(Q;e)"/ Inx, and then a — (1 — 5)/(l + In a). 
We also set wq = ln(l + In a). On the contour consisting of the lines ^{w) — wq + d, 
^(w) = ±25-'^/^, and 'R{w) = —2, one can show that |^| > |/|, and therefore f + g has only 
one root within the contour, by Rouche's theorem. Using Cauchy's theorem to express this 
root as an integral around the contour, we establish the convergence of (2a) by expanding 
the integrand as a series in f / g and integrating term by term [3]. For a < 1, the contour 
must remain the same as that for ct = 1. 
3. A New Expansion. — In view of the relation 

( X^l"^ 

$Q,(a;) = q;$i = aW 
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where W is the Lambert W function [6], we shall simplify our equations by considering 
only the case a = 1 from now on. By changing to the variable C = 1/(1 + cr), we obtain a 
new series for = $1 that converges on a wider domain than does (2a). 

Theorem 3. — With the preceding notation, W has the series development 
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(3a) Wix) = L,-L,+ J2'-j:i-ir--\^+-F + '^-'\ 
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and this is convergent for x >2. 

Proof. — Into (2c), we substitute a — l/( — 1 and obtain 

(36) T + e-'^^l + w-w/C = 0. 

To invert this using the Lagrange Inversion Theorem, we introduce the operator [tu^] to 
represent the coefficient of in a series expansion in w, and obtain 
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which can be rearranged to obtain the theorem. 

To prove convergence, we let f{w) = C(e~^ — 1 + w) + and g{w) = —w. On the 
rectangular contour bounded by the four lines = 2, '^{w) = ±2 and ^{w) = —1, it is 

simple to show that |/| < for all x G [2, e]. Hence the series converges there. Since (3a) 
is equivalent to (2a) for a; > e because of the relation 
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the theorem follows. 



4. Expansions using new variables. — Two further series developments can be 
obtained by introducing the variables Lr = ln(l — r) and rj = a/{l — r). 

Theorem 4. — With the preceding notation, W has the series development 
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Proof. — We set w — v — L^- m. (2c) and obtain, after rearranging, 
(46) 1 - e-^ + -^v = T^^r . 

1 — T 1 — T 

This equation has exactly the form of (2c) itself, and therefore the expansion for v can be 
obtained from {2d) by replacing a with a/{l — t) and r with aL^-f^l — r). The theorem 
then follows by rearrangement. 

The expansion (4a) converges more slowly than (2a), but when we transform it using 
the methods of theorem 3, we obtain a very rapidly convergent expansion, as we show in 
section 5. 

Theorem 5. — With the above notation, W{x) has the development 



(4c) w{x) = L,-L2-Lr + y ^L^v'^Yi-ir+^-'l^^'^ H \^ , 



Proof. — The proof follows exactly that of Theorem 3. 

The process of generating series in new variables can be continued. If w{a, r) satisfies 
(2c), then Theorem 4 is equivalent to the identity 

a a lii(l — t] 



(Ad) w{a, t) = -ln{l - t) + w 

which clearly can be applied repeatedly. 

5. Rate of convergence. — We consider the accuracy obtained by truncating each 
of the series (2a), (3a) and (4c) at N — 1 terms. Since the series are asymptotic series, 
the error terms for x large are respectively 0{L2 /L^) for (2a) and (3a) and 0{L2 /L\^) 
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for (4c), so (4c) is clearly better. In addition to being asymptotic, however, the series are 
absolutely convergent, and can be used for relatively small values of x. We observe that 
r = Lt- = at x = e, and hence the infinite sums in (2a), (3a) and (4c) are zero there. 
Thus any truncated series will be exact a,t x — e and asymptotically correct as x — > oo, 
implying that the error will have a maximum at some x > e. 

However, although (2a) is correct at a; = e, its derivative does not converge there. In 
contrast, (3a) and (4c) give finite sums at x = e for all derivatives. To put it another way, 
taking terms of (3a) or (4c) and expanding about x = e gives terms of the Taylor 
series for $i(a;) about x — e. Both (3a) and (4c) are much more accurate than (2a) near 
this point. Numerical experiments confirm these results. 

We conjecture that (2a) and (4c) converge for all a; > 1. 
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